Capillary condensation in disordered porous materials: hysteresis versus equilibrium 

behavior. 



O 

o 

(N 
C 

a 
oo 



E. Kierlik 1 , P. Monson 2 , M. L. Rosinberg 1 , L. Sarkisov 2 , and G. Tarjus 1 
1 Laboratoire de Physique Theorique des Liquidesu, Universite Pierre et Marie Curie, 4 Place Jussieu, 

75252 Paris Cedex 05, France 
2 Department of Chemical Engineering, University of Massachusetts, Amherst, MA 01003, USA 

We study the interplay between hysteresis and equilibrium behavior in capillary condensation of 
fluids in mesoporous disordered materials via a mean-field density functional theory of a disordered 
lattice-gas model. The approach reproduces all major features observed experimentally. We show 
that the simple van der Waals picture of metastability fails due to the appearance of a complex free- 
energy landscape with a large number of metastable states. In particular, hysteresis can occur both 
with and without an underlying equilibrium transition, thermodynamic consistency is not satisfied 
along the hysteresis loop, and out-of-equilibrium phase transitions are possible. 
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Capillary condensation of a gas inside a mesoporous 
material refers to the rapid change to a liquid-like state 
that occurs at a pressure (or chemical potential) lower 
than the bulk saturation value jj| . This phenomenon 
is often thought of as a shifted gas-liquid transition. The- 
oretical studies of fluids confined in single pores of ideal 
geometry have clarified the mechanism for such a shifted 
transition and introduced the concept of capillary criti- 
cality that describes the fact that the liquid- vapor critical 
point in a pore occurs at a lower temperature than in the 
bulk S [||. However, the situation in real mesoporous 
materials, such as porous glasses and silica gels, that 
consist of an interconnected network of winding pores of 
varying shape, curvature and size is not as clear. With 
the possible exception of fluids in aerogels of very high 
porosity (95 — 98%) for which a bona fide liquid-gas phase 
transition has been reported with no sign of hysteresis [Q , 
there is no direct evidence of a first-order phase transition 
characterized, say, by a jump in the adsorbed amount of 
fluid, nor of true criticality signaled by density fluctu- 
ations on very large length scales. The typical experi- 
mental observation is the presence of a hysteresis loop in 
the sorption isotherms, a loop that is reproducible and 
vanishes at a temperature lower than the critical tem- 
perature of the bulk fluid |Q j| U . The main questions 
raised by this phenomenology are: is there a true phase 
transition associated with capillary condensation in dis- 
ordered materials? What is the connection between this 
transition and the observed hysteresis? What are the 
sources of hysteretic behavior? No satisfactory answers 
have been provided so far. Indeed, in these systems that 
are locally heterogeneous, but macroscopically (statisti- 
cally) homogeneous, the statistical-mechanical methods 
that have proved very efficient in dealing with bulk fluids 
and fluids in simple external fields, namely computer sim- 
ulations, integral-equation theory and density functional 
theory, are plagued by many difficulties. 

In this letter, we present a mean-field density func- 



tional theory of a disordered lattice-gas model. This 
model, introduced in a previous work ||, incorporates the 
main physical ingredients characterizing fluids in disor- 
dered mesoporous media: preferential adsorption of one 
phase of the fluid, connectivity of the pore space, geomet- 
ric and energetic disorder, exclusion effect due to the ma- 
trix. (Note that the model could also be used to describe 
the demixing of binary mixtures in mesoporous materials 
Q.) The present approach allows us to study for the first 
time the interplay between out-of-equilibrium (hystere- 
sis) and equilibrium behavior associated with capillary 
condensation in a disordered matrix. We show that the 
theory captures the main aspects of the phenomenology 
of capillary condensation in real systems. The combi- 
nation of disorder and interconnectivity of the pore net- 
work drastically alters the picture of capillary conden- 
sation built upon the independent-pore model jlj Q 
or the simple van der Waals description of metastabil- 
ity because it generates a complex free-energy landscape 
characterized by a large number of metastable states. In 
particular, we show that (i) hysteresis can occur with or 
without an underlying equilibrium phase transition, (ii) 
the disappearance of the hysteresis loop is not associated 
with capillary criticality, (iii) thermodynamic consistency 
is not satisfied along the hysteresis loop, and (iv) collec- 
tive out-of-equilibrium phenomena such as macroscopic 
"avalanches" are possible. 
The model is described by the following Hamiltonian, 
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where r, = 0, 1 and (I — ?/;) =0,1 denote the fluid and 
matrix occupancy variables, respectively, and the sums 
run over distinct pairs of nearest-neighbor (n.n.) sites. 
Although different disordered microstructures of the ma- 
trix could be considered pj , we choose here the simplest 
nontrivial random matrix. The model is thus 
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specified by two parameters, the average matrix density 
p m that fixes the porosity (equal to 1 — p rn ) and the 
ratio of the matrix-fluid over the fluid-fluid interactions, 
y = w m f /w f f , that determines the "wettability" . In the 
following we only consider attractive matrix-fluid inter- 
actions H and we set p m = 0.25. We use a 3D bcc lattice 
with linear size L = 8, 12, 16, 24, 48 for which we consider 
many different configurations of the random matrix, from 
1400 for L = 8 to 5 for L = 48. Most of the results are 
illustrated here for L = 48 and typical matrix configura- 
tions, but when necessary, e.g., to check the existence of 
phase transitions, we have performed a finite-size analysis 
after averaging over all matrix configurations. 

The mean-field density functional theory (or equiva- 
lently for a lattice, the local mean-field theory) starts 
with the formulation of a free-energy functional of the 
fluid density field for a given matrix realization {rji}: 

f apt}) = 4 lnpi + ~ p*) ln (^ ~ p^ 

i 

- w ff Y p ip i ~ Wm f Y t^ 1 ~ + p^ 1 ~ ' ( 2 ) 

<ij> <ij> 

where = 1/(/sbT) and pi —< 7ft Tj > is the average 
fluid density at site i. (The overall fluid density is then 
Pf = (l/N)J2i Pi, where N — 2L 3 is the total number 
of sites.) For a given chemical potential p, minimization 

of fl({pi}) = F({pi}) - fJ-J2iPi w - r - t - tne Pi' s Pro- 
vides the grand potential fl of the adsorbed fluid. The 
corresponding equations for the fluid density on each site 
are 

Vi 



1 + exp 
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P + Ej/i ( w ffPj + w mf (! - Vj)) 



(3) 



where the sum is over all n.n. of site i. The above set 
of nonlinear coupled equations has been solved by means 
of a simple iteration algorithm according to two differ- 
ent protocols: (1) to mimic the experimental procedure, 
we follow continuously the solutions under small varia- 
tions of the chemical potential p (typically, Ap = 10~ 3 ) 
and (2) to search more exhaustively the solutions of Eq. 
(3) for a given p, we repeat the iteration procedure with 
a large number of initial conditions (typically, 10 2 ) cor- 
responding to uniform fillings of the lattice at different 
overall fluid densities (i.e., p^ = p^)- 

We show in Fig. 1 typical isotherms obtained by in- 
creasing continuously p from — oo (adsorption branch) 
and decreasing continuously p from +oo (desorption 
branch). They look qualitatively similar to the exper- 
imental isotherms of fluids adsorbed in porous glasses 
and silica gels [0 |@| @, with an asymmetric hystere- 
sis loop (corresponding to the H2 shape in the IUPAC 
classification) characterized by a steep desortion branch 



CX 0.4 - 




8/.i 



FIG. 1. Theoretical sorption isotherms for y = 1.5 at 
T = 0.6,0.8,0.95,1.1 and 1.25. (The mean-field bulk criti- 
cal point is at T = 2.) 



and a smoothly increasing adsorption branch. The so- 
called "hysteresis phase diagrams" that are constructed 
from the points at which the loop opens and closes are 
shifted to lower T and higher pf compared to the bulk- 
fluid coexistence curve, and their (asymmetric) shape is 
also quite comparable to what is obtained from experi- 
mental data @] § §. 

An important test for the relevance of any model of 
capillary condensation in disordered porous materials is 
the shape of the so-called "scanning curves", i. e., curves 
that are obtained by reversing the sign of the evolution 
of p at different points along the adsorption or the des- 
orption branches M. It has been stressed that these 
curves are not properly reproduced by the widely used 
independent-pore model |l| [|| • An illustration of the 
scanning curves on adsorption and desorption obtained 
within the present theory is shown in Fig. 2. They look 
strikingly similar to those observed experimentally, with 
an upward curvature on adsorption and a downward cur- 
vature on desorption. When varying the chemical poten- 
tial up and down along different paths we find a hierarchy 
of inner scanning curves that display the "return-point 
memory" effect, in close analogy to the hysteretic behav- 
ior displayed by a variety of systems [|l0| [p] . 

The above results confirm that the present approach 
reproduces qualitatively the phenomenology associated 
with capillary condensation in porous glasses and silica 
gels. This indicates that, despite its simplicity, the model 
does capture the main physical features of fluids adsorbed 
in real mesoporous systems and that equilibration via 
thermally activated processes, processes that are absent 
in the mean-field description, does not substantially mod- 
ify the picture. Indeed, the implicit assumption behind 
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FIG. 2. Theoretical desorption (a) and adsorption (b) scan- 
ning curves for y = 1.5 and T — 0.8. 

the direct comparison of the theoretical sorption curves 
with the experimental ones is that the system does not 
have time to equilibrate at constant T and p and that ad- 
sorption/desorption only proceeds under the influence of 
chemical-potential changes. (This assumption may break 
down when the perturbation induced by the matrix is 
small as in very dilute aerogels.) 

Hysteresis and scanning curves are manifestations of 
out-of-equilibrium behavior. This is well understood 
for instance in systems like athermal martensites and 
low-T ferromagnetic materials whose transformations are 
driven by an external field; these systems are well de- 
scribed by the T = limit of models with quenched dis- 
order |Kj| (models somewhat related to the lattice-gas 
studied here). On the other hand, in the case of cap- 
illary condensation one often assumes that hysteresis is 
also the signature of a true equilibrium phase transition 
with two possible metastable states, gas and liquid, as 
in the standard van der Waals loop for bulk fluids and 
fluids confined in a single pore. However, this simple van 
der Waals picture fails for the present system because 
this latter has (in the region where capillary condensa- 
tion occurs) a complex free-energy landscape with a large 
number of metastable states, as can be seen from study- 
ing the solutions of the local mean-field equations, Eq. 
(3). For each T, p and each matrix realization, we have 
investigated 50 to 100 different initial conditions to Eq. 
(3). When plotted as in Fig. 3 on the pf — p diagram 
(for a typical matrix configuration of the L — 48 system) , 
all solutions fall inside the major hysteresis loop that is 
found to coincide with the curves of extremal solutions 
obtained from an initially empty lattice (lower branch) 
and an initially filled lattice (upper branch). We suspect 
that this picture corresponds to the situation encountered 
in real disordered porous materials (a similar behavior 



FIG. 3. Multiplicity of local mean-field solutions obtained 
from solving Eq. (3) with many different initial conditions for 
each value of jj,: (a) y — 1.5 and T = 0.8, (b) y = 1 and T = 1. 
The solid lines represent the equilibrium curves obtained by 
connecting the states of lowest grand potential. 

has been observed for other disordered microstructures 
1). 

The most significant finding of our study is that hys- 
teresis associated with capillary condensation in disor- 
dered porous media occurs with as well as without an 
underlying equilibrium phase transition, the existence of 
this latter depending on the strength of the perturbation 
induced by the matrix. We have computed the equilib- 
rium properties of the model by means of the weighted 
mean- field approach at each p and T, densities of 
extensive quantities such as the fluid density pj are cal- 
culated as Boltzmann weighted sums over the different 
solutions of the mean-field equations (i.e., each solution 
a is weighted by e _/3S1 / ^2 a e~^ n ). The average over 
the various matrix realizations is taken afterwards. For 
a large matrix-fluid interaction, e.g., for y = 1.5, there 
is no equilibrium phase transition: this is illustrated in 
Fig. 3a where the equilibrium isotherm at T = 0.8 is 
perfectly smooth (it stays so at lower temperatures) . On 
the other hand, for y = 0.5, the system undergoes a true 
thermodynamic transition at low T as it becomes iso- 
morphic to the site-diluted Ising model ^ (p m = 0.25 
is above the percolation threshold). For the intermedi- 
ate value y = 1, equilibrium isotherms may display steep 
portions, as illustrated in Fig. 3b for T = 1 (however, 
a finite-size scaling analysis of the results for L — 8,12 
and 16 seems to indicate that this isotherm is still super- 
critical). The capillary critical point certainly occurs at 
a T significantly below the temperature Th ys t at which 
hysteresis first appears. This result is at odds with the 
behavior predicted for a fluid confined in a single pore j|] 
H and is a consequence of the quenched disorder; on the 
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other hand, the fact that true capillary criticality may 
still occur in disordered porous media is at odds with the 
prediction based on the independent-pore model |D @] 
and results from the connectivity of the void space 
accessible to the adsorbed fluid. 




Q. 0.4 




-5.5 -5 



/i./w [r 



FIG. 4. Check of thermodynamic consistency 
((dfl/dfi) T = —Npf) along the adsorption (a), desorption 
(b), and equilibrium (c) isotherms (y — 1.5, T = 0.8). Filled 
symbols: average fluid density obtained from the mean-field 
solutions. Open symbols: quantity obtained by differentiating 
the corresponding grand potentials. 



The breakdown of the simple van der Waals picture of 
metastability in the present problem shows up distinctly 
when considering thermodynamic consistency along the 
sorption isotherms. Since the theory provides both the 
grand potential and the fluid density, one can study the 
validity of the Gibbs adsorption equation, (dQ/d[i) T = 
—Npf, along particular isotherms. As shown in Fig. 4 
a,b, the equation is not satisfied, neither along the ad- 
sorption isotherm nor along the desorption one. This 
results from the fact that the system often jumps from a 
given grand potential minimum to another along these 
two isotherms and then looses thermodynamic consis- 
tency. An important consequence is that the thermo- 
dynamic integration procedure used for building equilib- 
rium phase diagrams is no longer valid for this model. 
For instance, in the case y = 1.5, this latter procedure 
predicts a capillary phase diagram whereas, as already 
discussed, no equilibrium phase transition takes place. 
We believe this may also be true more generally for other 
models of fluids in disordered porous materials. Note 
that the Gibbs adsorption equation is obeyed along the 
equilibrium isotherm, as illustrated in Fig. 4c. (Interest- 



ingly, in the cases studied here, the isotherm calculated 
from the weighted mean-field approach is virtually indis- 
tinguishable from that obtained by taking at each \x the 
solution with the lowest grand potential.) 

As a final remark, we mention that the mean-field den- 
sity functional theory also predicts the occurrence of out- 
of-equilibrium phase transitions. Similarly to the macro- 
scopic avalanches observed in low-T ferromagnetic mate- 
rials [|l0| , the adsorption and desorption branches of the 
major hysteresis loop may display sharp jumps. This is il- 
lustrated in Fig. 1 where one can see that the desorption 
isotherm has a discontinuity at the lowest temperature 
shown. Due to the asymmetry of the random chemical 
potentials exerted by the matrix on the fluid, the jump 
always occurs at a higher T in the desorption branch 
than in the adsorption branch. Whether or not this out- 
of-equilibrium transition, with its associated critical be- 
havior |l(]] , can be observed in capillary condensation of 
fluids in real mesoporous materials depends on the effi- 
ciency of the thermally activated processes. 
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